Collective modes in quantum electron glasses and electron-assisted hopping 
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We study electronic transport in Anderson insulators with strong Coulomb interactions in dimen- 
sions d > 2. Close to the metal insulator transition where the single particle localization length 
is much larger than interparticle-distance, the interactions lead to a strongly correlated quantum 
glass phase. Even though single particle excitations are localized and the system is insulating, there 
are collective electronic modes which remain delocalized down to parametrically small energies. 
These collective excitations serve as a continuous bath which can provide the activation energy for 
variable range hopping transport. This circumvents the energy conservation problem arising when 
only discrete particle-hole excitations are present. In contrast to the weak and material-dependent 
phonon-assisted hopping mechanism, the activation by an electronic bath leads to a nearly universal 
prefactor e'^/h of the Efros-Shklovskii conductance, as is observed in many recent experiments. 

PACS numbers: 71.23.An, 71.30.-(-h, 72.20.Ee, 64.70.pc, 81.05.Gc 



Close to the metal-insulator transition (MIT), 
Coulomb interactions turn Anderson insulators into 
strongly correlated and frustrated quantum systems. 
Their transport properties have been a puzzle for a long 
time and are still lively debated. On the experimen- 
tal side there is a longstanding discrepancy between the 
observed, almost universal pre-exponential factor of the 
hopping conductance [ij (of order e^/h in two dimen- 
sions), and the much smaller and material-dependent 
value predicted by the standard theory based on the acti- 
vation of electrons by phonons [3| ■ The universality of the 
prefactor lead to speculations about a purely electronic 
activation mechanism 0, i4| , but a consistent microscopic 
theory remains elusive. In fact, the situation became 
even less clear recently, when several authors showed that 
purely electron-assisted transport is impossible because 
of localization in Fock space, if interactions are weak and 
sufficiently short ranged The paper Q claims, how- 

ever, that this dichotomy is resolved by the long range 
nature of Coulomb interactions, rendering Refs. [1, @] in- 
applicable to physical systems in d > 2, but it leaves the 
question of the prefactor open. So does the theory of 
electron-assisted variable range hopping constructed for 
disordered Luttinger liquids in d = 1 [8|, which, in ad- 
dition, relies on the presence of a thermal phonon bath 
and thus is in fact not purely electronic in nature. 

However, even without invoking long-range interac- 
tions, many body localization is unlikely to occur if the 
interactions are strong enough, as suggested by exact di- 
agonalization studies in one- and two-dimensional mod- 
els H, 0. This is indeed what should be expected as 
the MIT is approached and single electron wavefunc- 
tions start to overlap: The strength of Coulomb inter- 
actions decreases less rapidly than the level spacing in 
a localization volume. This drives the electronic sys- 



tem to a strong coupling regime, where perturbation 
theory breaks down and a collective, quantum glassy 
state emerges. In this paper, we analyze this regime and 
show that, even though the system remains insulating, 
it possesses a continuum of delocalized collective elec- 
tronic modes whose coupling to hopping electrons entails 
a nearly universal pre-exponential factor in the hopping 
conductivity, of the order observed in experiments. 

We consider interacting electrons in d > 2, with dis- 
order such that the single-particle localization length ^ 
is much bigger than the average interparticle spacing, 
i.e., n^'' ^ 1, where n is the electron density. As a ba- 
sis for the interacting problem we choose the localized 
eigenfunctions tpi (r) of the selfconsistent disorder poten- 
tial, which accounts for the screening due to the "frozen" 
electron density common to all mctastable states of the 
system. In this basis the Hamiltonian reads 

H ^ ^'^UiJijUj + '^{ei- |J.)n^ (1) 

i^j i 

where ^ is the chemical potential and Srii = rii — {{n-i)) 
is the deviation of the occupation number rii from the 
"frozen" background occupation ((n;)). The denote 
the energies of the non-interacting problem character- 
ized by a bare density of states vq. The typical level 
spacing in a localization volume is thus = l/{£,'^i^o). 
Jij and t^.j , V^^'' are the diagonal and off-diagonal matrix 
elements of the Coulomb interactions, respectively. Due 
to the random charge density pattern \tpf{r)\ associated 
with the level i, the Coulomb interaction between two 
states within the same localization volume is a random 
quantity of order \Jij \ ^ J = / k^, after subtraction of 
an irrelevant average repulsion. Here, k is the dielectric 
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constant of the host. The strength of Coulomb interac- 
tions is measured by the ratio 



(2) 



z = 

which diverges as the MIT is approached — * oo). Once 
the wavefunctions start overlapping, z rapidly becomes 
large, z ^ I, leading to a highly frustrated glassy prob- 
lem. This situation is opposite to the limit z <^ 1 ana- 
lyzed in Refs. which allowed for a perturbative treat- 
ment starting from the non-interacting ground state and 
lead to the prediction of many body localization. 

The first line in ^ defines a classical problem equiv- 
alent to an Ising spin glass with random fields. Note 
that z^ 1 plays the role of a large effective coordi- 
nation number [lol . [ll[ , which allows for controlled ap- 
proximations. The classical Hamiltonian possesses many 
mctastable states, i.e., low energy configurations rii which 
are stable with respect to the rearrangement of few elec- 
trons. The Hartree interactions with other electrons 
strongly modify the distribution of the effective energies 
of the levels, = — /i + Jij'^'j-, creating a lin- 
ear pseudogap similar to the Efros-Shklovskii Coulomb 
gap l3, la| which extends up to the large energy scale 
Ec ~ zJ The Coulomb interactions reshuffie the 

single particle energies, and in a given localization vol- 
ume the average level spacing at the Fermi level is renor- 
malized to J ^ 5^. At temperatures of order J one thus 
expects the onset of activated hopping conductivity. 

In the presence of quantum ffuctuations due to the 
off-diagonal terms of the Coulomb interactions, stability 
arguments show that the pseudogap in the quasiparticle 
energy distribution remains intact. Further, according to 



standard estimates [15| the matrix elements t'^.j and 
decrease with bare energy transfer on average like 
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J$((|6,-6fc|)/J), (3) 

J<I>((|6,-e,|)/<5)$((|6,~£n)/'5), (4) 



with - x"^/"^ for x < {i/(.Y\ where £ < ^ is the 

mean free path. 

The large effective coordination number and the 
slow decay of matrix elements ensure that a given level 
interacts with many neighbors. This allows us to neglect 
fermion loops and adopt a bosonic pseudo-spin 1/2 de- 
scription rii ^ si + 1/2, (cj, Cj) ^t'~ ^ ^ii*^ to look for 
a self-consistent mean-field description of quantum fiuc- 
tuations with the effective action (similarly to Refs. [l6l |) 



5, 



cff 



dr 



/3 



-f^/ \ dra+{r')G,{r'-T)<j-{r). (5) 
j Jo Jo 

The kernels G,; (r) are due to the virtual hopping of elec- 
trons to other sites and back. Thus, they describe in- 
ertial; non-dissipative dynamics, which does not lead to 



significant level broadening, as we will verify below. In 
order to gain insight into the collective modes of this 
system we proceed similarly as in the TAP approach to 
classical spin glasses ll|, |l3[, which is a good approxi- 
mation for z ^ 1. The effective field acting on spin i 



IS 



hj Xj J 



(6) 



where the "magnetizations" = (erf) are averages re- 
stricted to a given metastable state. The last term in 
(l6|) describes the Onsager reaction of the environment 
of neighboring spins, where Xj is the susceptibility of 
spin J. An exact treatment would require to formu- 
late self-consistent dynamic TAP equations involving the 
dynamic susceptibilities Xjk^ ~ ''"') [l3|- However, the 
short-time dynamics describing collective modes can be 
captured by the much simpler static approximation [l8l | 
([6|); which can be derived as saddle point equations 
dJ-'{{mj})/dmi = of an approximate static energy func- 
tional following Ref. [lit . 

The equations ^ are closed when supplemented by a 
relation determining the magnetization m.i in the pres- 
ence of a field hi. This can be achieved by realizing that 
quantum fluctuations act like effective transverse fields 
of order icff ^ Ec for the spins. We have estimated toff 
by calculating the magnetizations of sites in strong lo- 
cal fields with a perturbative calculation based on the 
Hamiltonian ([T]), and comparing them with the corre- 
sponding expression resulting from an effective transverse 
field. The reduction in magnetization is due to hybridiza- 
tions with other levels, predominantly mediated by the 
electron-electron scattering term y^^'. Since variations 
of the fluctuation strength are inessential for the physics 
we approximate the effective transverse field for simplic- 
ity by a constant magnitude toff Ec. This provides 
us with the sought relation hi = 2tcsrai/[l/A — mf]^/^ 
between field and magnetization, and Xi — — 

Yj^ [l/4— mf]^^^ for the susceptibility. 

The large magnitude of quantum fluctuations is due to 
the marginal decay of (|3I4[) . This ensures a large number 
■N'act ~ ^tsl^\ ~ -^^ electrons per localization volume 
which participate significantly in the quantum dynamics 
of the glass, justifying our mean-field approximations. 

In order to assess our approximate TAP analysis, we 
have tested it on the case of a fully connected spin model 
with random Gaussian couplings Jij {Jf^ = J'^/N) in a 
transverse field t, but without random fields (e,; = 0). 
One finds that the instability of the paramagnetic TAP 
solution, rrii = 0, occurs precisely at the same critical 
transverse field tc = 2y/N J as predicted by the static ap- 
proximation of the replica approach [20| . which suggests 
that the two approximations are equivalent. Both under- 
estimate the strength of fluctuations by approximating 
the Onsager response by the full adiabatic susceptibility 
of spins instead of their smaller dynamic response. 
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Solutions {rrii} of Eqs. ([6]) represent metastable states 
of the quantum glass. Collective excitations around these 
local minima are governed by their potential energy en- 
vironment which is captured by the Hessian 



Hi 



dnij 



dmidnij 

where we have defined the 'locator' L 



LiSij , 



(7) 



X, 



Similarly as in classical glasses |2l| and in mean field 
quantum glasses [l3, O'^^ expects that for J ^ toff, 
the spectrum of Ti is gapless, with a density of eigenvalues 



P-h(A) 



TT 
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-1/2 



J2 



(8) 



for A < tcS- The gaplessness within the glass phase re- 
flects the marginal stability of the quantum glass which 
self-organizes into states {mi\ characterized by a critical 
Onsager response, '^ijXj — 1; ^"^^ abundance of 
collective low energy modes. 

For the purpose of understanding the low energy spec- 
trum, the Hessian ([7]) can be restricted to the "quantum 
active" spins with < ieS, recalling that any spin i in- 



teracts substantially via Jy- with A/'a, 



active spins. 



The eigenvectors of random matrices with this structure 
have been analyzed in One finds that there is a 

parametrically small mobility edge Ac ^ Jz^^l"^ in the 
spectrum, above which all modes are delocalized. 

The quantum dynamics of these collective excitations 
is non-dissipative at low energies, as we will check below. 
We thus assume them to behave like harmonic oscilla- 
tors with restoring force A, and a mass M ~ inde- 
pendent of the mode (in the regime of interest A < toff), 
as will be justified later. The collective modes are thus 
characterized by a frequency ujx = (Atoff)^^^ and mean 
square displacement {x^) ^ ~ tcff/lfj^l- From the density 
of restoring forces ([5]) we then easily obtain the spectral 
function of collective modes (for u) <^ toff) 



X"{U)^{X^)^ d\p{\)5{LU-L,{\)) 



(9) 



which, remarkably, turns out to be independent of toff . 
Eq. ^ is equivalent to a local dynamic density correlator 
with long time asymptotics ^^'^ {p{r,T)p{r,0)) ^ l/(Jr)^. 
We have adjusted the oscillator mass M in such a way 
that this correlator scales like A/^ct for r ^ AI, reflecting 
that all active spins take part in the dynamics at these 
short time scales. The power law decay of density 
correlators and linear spectral functions x"{^) \^\ ai'e 
well-known features of quantum glasses with undamped 
dynamics at and below the glass transition Here 
we obtain this result with the physically intuitive inter- 
pretation of collective modes as underdampcd harmonic 
oscillators governed by the semicircular distribution ([5]). 



We have verified that neither particle-hole excitations nor 
the coupling to other collective modes provide substantial 
friction, so that the description in terms of underdamped 
oscillators is indeed correct. 

We emphasize that the collective modes with energies 
above wioc = ujx^ ^ Jz^^l^ are delocalized, and hence 
form a bath with continuous spectrum. In linear approx- 
imation low energy modes with lo < wioc are localized, 
but it is possible that cubic couplings among them in- 
duce delocalization which might lower the edge ujc where 
the spectrum of collective modes turns discrete. Even 
though it seems unlikely, we cannot rule out that tUc ^ 0. 

The collective modes are globally neutral charge den- 
sity oscillations which cannot contribute directly to the 
electrical conductivity. However, they can activate sin- 
gle particle transitions close to the Fermi level, leading 
to eZectron-assisted hopping. Note that the continuity of 
the spectrum is crucial to ensure energy conservation in 
an elementary hopping process. In contrast, particle-hole 
excitations are discrete and thus ineffective as a bath for 
inelastic hopping processes 0, 0| • 

It is not immediately obvious that in the presence of 
quantum fluctuations and collective low energy degrees 
of freedom the system remains an insulator. However, 
we have checked that single particle-excitations remain 
localized. Further, at low temperature the level broad- 
ening is smaller than the level spacing at the Fermi level 
in a localization volume. Indeed, decay rates of low-lying 
states are very small since they involve transitions with 
strongly suppressed matrix elements. On the other hand, 
pure dephasing due to the coupling of single particle lev- 
els to the collective bath leads to a level broadening of 
the order of I/T2 ^ T. This is negligible at tempera- 
tures T <ti J where Efros-Shklovskii type variable range 
hopping sets in We also note that the large pa- 

rameter z allows us to neglect the small overlap a single 
particle excitation has with a collective mode. 



Following standard arguments of hopping theory [24| , 
and replacing phonons by the electron bath with spectral 
function ([9]), one obtains a percolation problem of Miller- 
Abrahams resistors with elementary conductances 

""^"^nF(e^)[l - nF(£j)],(10) 



ff(A£y)e" 



where is the Fermi distribution and Ae^ 



is the energy to transfer an electron from site i to j . For 
electron-assisted activation we find g{Ae) ^ XTS^' 
which depends only on the strength of Coulomb in- 
teractions, J. This is in marked contrast to phonon- 
assisted hopping which leads to a much smaller prefac- 
tor g proportional to the material-dependent electron- 
phonon coupling constant. The bulk conductivity fol- 
lows from percolation arguments for the network of el- 
ementary resistors with conductances ([T0| [23 |. For d- 
dimensional systems this yields the prediction 



ff(£^hop)r"' (e/i?hop 



-iJ/T) 
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where the energy and distance of typical bottleneck hops 
are Ehop/T ~ i?hop/C ~ {J/TY^^. With the percolation 
exponents I'd « 1.34 (0.9) for d = 2 (3), respectively, one 
finds a weak T dependence of the prefactor, with yu^ = 
3/2 -(l + iyd)/2 = 0.33 (0.55). 

Note that for d = 2, our theory predicts a nearly uni- 
versal prefactor of order e'^ /h (as obtained by extrapo- 
lating conductivity plots to J/T 0), while for d = 3 
it is of the order of the minimal metallic conductivity. 
Such behavior has been reported in many recent exper- 
iments the prefactor being consistent with small or 
vanishing > 0. Note however, that the latter is sensi- 
tive to assumptions about the matrix elements involved 
in elementary hopping processes, and is therefore less ro- 
bust than the near universality of the prefactor. 

We emphasize that the above derivation is valid only 
for temperatures such that the activation energy of bot- 
tleneck resistors exceeds the " mobility edge" in the spec- 
trum of collective modes, i?hop(2^) > (^c- At lower tem- 
peratures the conductivity becomes simply activated of 
the form a cx exp(— Wc/T), but still with a prefactor of 
electronic origin. At even lower temperature, the phonon 
bath will take over and restore variable range hopping, 
however with a significantly smaller prefactor. Apart 
from the clear signature in the conductivity, the pres- 
ence of delocalized collective electronic modes is expected 
to show experimentally in their contribution to thermal 
conductivity [i^l, which should decrease substantially as 
the temperature drops below ujc- 

In the derivation of the properties of the quantum glass 
we have tacitly assumed that the quantum fluctuations, 
as measured by tes, are smaller than the critical value 
tc ~ Ec which would melt the glass [11] and lead to a 
metallic state with delocalized quasiparticle excitations. 
While it is difficult to imagine that this could result from 
repulsive Coulomb interactions, we cannot directly rule 
out this possibility at very large z. However, we can in- 
fer the existence of a stable quantum glassy regime with 
overlapping electron wavefunctions (z > 1), from the es- 
tablished existence of a glass phase [2^ in the nearly 
classical regime of well localized electrons in moderate 
disorder where z ^ I. This conjecture is borne out 
by experimental observation of quantum glassy behav- 
ior in amorphous insulators with high carrier density 27 1 . 
Those systems can be taken into a regime of large z where 
glassy effects in fact become stronger [2^ . The puzzling 
density and temperature dependence of glassy relaxation 
reported in [2^ might be related to the emergence of low 
energy excitations discussed in this paper. 
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